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Abstract. We present a new search for the Ks — > Sir decay performed with the KLOE 
detector operating at the DA$NE <j!>-factory. The Ks mesons were tagged via registration of 
Kl mesons which crossed the drift chamber without decaying and interacted with the KLOE 
electromagnetic calorimeter. The Ks — > 3ir° decay was then searched requiring six prompt 
photons. To suppress background, originating from fake Ks tags and Ks — ► 2tt° decays 
with additional two spurious clusters, we have performed a discriminant analysis based on 
kinematical fit, testing of the signal and background hypotheses and exploiting of the differences 
in kinematics of the Ks decays into 2tt° and Sir . In a sample of about 1.7T0 9 <f> — > KsKl 
events we have found no candidates in data and simulated background samples. Normalizing 
to the number of Ks — > 2n° events in the same sample, we have set the upper limit on the 
K s -> 3tt° branching ratio BR(K S -> 3tt°) < 2.6 • 1CT 8 at 90% C.L., five times lower than 
the previous limit. This upper limit can be translated into a limit on the modulus of the 77000 
parameter amounting to 77000 1 < 0.0088 at 90% C.L., improving by a factor two the latest direct 
measurement. 



1. Introduction 

Since the first discovery of the CP-violating neutral kaon decay in 1964, there has been made 
a big effort to describe the CV symmetry breaking within the Standard Model. The favoured 
theoretical framework was provided in 1973 by Kobayashi and Maskawa, who pointed out that 
CV violation would follow automatically if there were at least six quark flavours. At present the 
main experimental effort is focused on the neutral B and D meson system studies. However, 
there are still several interesting open issuses in the kaon physics, which can contribute to our 
better understanding of the CV violation mechanism [1]. For example, Kg decays to |7r + 7r~7r ) 
and |37T°), requiring CV violation, are still poorly known. The |7r + 7r~7r°) final state can be 
produced in neutral kaon decays with isospin I = 0,1,2, or 3. The 1 = and 1 = 2 states 
have CV = 1, and Ks can decay into them without violation of the CV symmetry. However, 
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they are expected to be strongly suppressed by centrifugal barrier effects [2]. For the 1 = 1 
and 1 = 3 states there is no centrifugal barrier and CV = —1, so K$ decay requires violation 
of this symmetry. Anyhow the two kinds of final states can be separated by the analysis of 
the 7r + 7r~7r Dalitz plot. This allows for determination of the Kg to Kl decay amplitude ratio 
?? + _o = a{k S l ^Z+1-1°] - € + e +-o> where e indicates the K s CV impurity and e' + _ is the 
contribution of the direct CP-violating term. In the case of \ir ir°ir ) final state, only isospin 
I = 1 or I = 3 is allowed, for which CV = — 1. Therefore, the K$ — > 3tt° decay is a purely CV 
violating process, for which one defines the analogous ratio: 77000 = ^(^Z^o^o^o) — e + £ooo- 

The present knowledge about r/_| an d ??ooo is poor mainly due to very low decay rates for 

the Ks — > 3-7T decays. The current value of the Ks — > ir + ir~ir° branching ratio amounts to 
BR(K S -> tt+tt-tt ) = ( 3.5+^ ) • 10~ 7 [2j, and the K s -> 3vr° has been never observed. The 
best upper limit on this decay branching ratio amounts to BR(Ks —5- 3ir°) < 1.2- 10~ 7 [3j, while 
the prediction based on Standard Model is equal to about 2 • 10 -9 [3j. 

In this article we briefly describe the search of the Ks — > 3n° decay based on the 1.7 fb _1 of 
integrated luminosity gathered with the KLOE detector operating at the ^-factory DA<I>NE of 
the Frascati Laboratory [5]. 

2. The KLOE detector 

The KLOE experiment operated from 2000 to 2006 at DA^NE, the e + e~ collider optimized to 
work with a center of mass energy around the 4> meson mass peak: ^/s = 1019.45 MeV [B]. Equal 
energy positron and electron beams collided at an angle of tt — 25 mrad producing eft mesons 
nearly at rest, which then decayed mainly to K + K~ (49%) and KsKl (34%) final states. 
The decay products were registered using the KLOE detection setup presented schematically 
in Fig. [TJ The detector consists of a large cylindrical Drift Chamber [7] surrounded by a lead 
scintillating fiber Electromagnetic Calorimeter [8], both immersed in an axial 0.52 T magnetic 
field produced by a superconducting solenoid. The beam pipe at the interaction region is a 
sphere with 10 cm of radius made of a Beryllium-Aluminum alloy of 0.5 mm thickness. The 
drift chamber, 4 m in diameter and 3.3 m long, has 12582 all stereo drift cells with tungsten sense 



wires. To minimize the Kl regeneration, multiple Coulomb scattering and photon absorption 
the chamber is constructed out of carbon fiber composite with low-Z and low density, and uses 
the gas mixture of helium (90%) and isobutane (10%). It provides three-dimensional tracking 
with resolution in the bending plane about 200 fim, resolution on the z-coordinate measurement 
of about 2 mm and of 1 mm on the decay vertex position. Momentum of a particle is determined 
from the curvature of its trajectory in the magnetic field with a fractional accuracy <J p /p = 0.4% 
for polar angles larger than 45° The KLOE electromagnetic calorimeter covers 98% of the 
solid angle and is composed by a barrel and two endcaps, for a total of 88 modules. Each 
module is built out of 1 mm scintillating fibers grouped in cells of 4.4x4.4 cm 2 and embedded 
in 0.5 mm lead foils, and it is read out from both sides by set of photomultipliers. The particle 
energy deposits are obtained from the photomultipliers signal amplitude, while the arrival times 
and particles impact points are determined from the spatial coordinates of the fired calorimeter 
cell and the difference between times measured at both ends of the calorimeter module. The 
accuracies of energy and time determination are parametrized as cje/E = 5.7%/ \JE (GeV) and 
at = 57 ps/ y/E (GeV) ® 100 ps, respectively. The trigger [9] uses both calorimeter and chamber 
information. In the search for the K$ — > 3ir° we have used, however, only the calorimeter 
information requiring two energy deposits with E > 50 MeV for the barrel and E > 150 MeV 
for the endcaps. 

3. Data analysis 

At KLOE kaons arising from the (ft decay move at low speed with their relative angle close to 
180°. Therefore, observation of a Kl (K$ ) decay ensures the presence of a K$ {Kl ) meson 
travelling in the opposite direction. About 50% of Kls reach the calorimeter before decaying, 
which provides a very clean Kg tag by the Kl interaction in the calorimeter (if^— crash). It 
is identified by a cluster with polar angle 40° < cr < 140°, not associated to any track, with 
energy E cr > 100 MeV and with a time corresponding to a Kl velocity (3* ~ 0.2 in the (ft rest 
frame. The average value of the e + e~ center of mass energy y/s is obtained with a precision of 20 
keV for each 200 nb _1 running period using large angle Bhabha scattering events [BJ- The value 
of yfs and the -crash cluster position allows us to obtain, for each event, the direction of the 
Kg with an angular resolution of 1° and a momentum resolution of about 2 MeV. The search for 
the Ks — > 37T° —7- 67 decay is then carried out by the selection of events with six photons which 
momenta are reconstructed using time and energy measured by the electromagnetic calorimeter. 
Photons from the Ks decay are reconstructed as neutral particles that travel with a velocity 
(3=1 from the interaction point to the KLOE calorimeter (prompt photons) . In order to retain 
a large control sample for the background while preserving high efficiency for the signal, we 
keep all photons satisfying E 1 > 7 MeV and | cos6*| < 0.915. Each cluster is also required to 
satisfy the condition |£y — FL,/c\ < min(3.5o"t,2 ns), where t 7 is the photon flight time and R 
the path length. The photon detection efficiency of the calorimeter amounts to about 90% for 
E 1 = 20 MeV, and reaches 100% above 70 MeV. To determine the Ks — > 3tt° branching ratio 
we have counted also the Ks — > 2-7T events which are selected requiring 4 prompt photons, and 
have served as normalization sample. For both channels the expected background as well as the 
detector acceptance and the analysis efficiency are estimated using the Monte Carlo simulation 
tool of the experiment |10| . All the processes contributing to the background were simulated 
with statistics twice larger than the data sample. Moreover, for the acceptance and the analysis 
efficiency evaluation a dedicated Ks 37r° signal simulation was performed. After Ks tagging 
we have counted 76689 events with six prompt photons. For these events a further discriminant 
analysis was performed to increase the signal to background ratio. 

The first analysis step aimed to reject fake Ks tags. The distributions of reconstructed Kl 
energy E cr and velocity (3* for the selected data sample and background simulations are shown 
in Fig. O In the (3* distribution, the peak around 0.215 corresponds to genuine Kl interactions 
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Figure 2. Distributions of the reconstructed Kl energy (E cr ) and velocity in the (ft center of 
mass frame for all events in the six photon sample. Black points represent data, while the 
MC background simulation is shown as red histogram. 

in the calorimeter, while the flat distribution mainly originates from (ft ~^ K$Kl — > {K$ — > 
tt + it~,Kl — > 37T°) background events. In this case one of the low momentum charged pions 
spirals in the forward direction and interacts in the low-/? insertion quadrupoles. This interaction 
produces neutral particles which simulate the signal of Kl interaction in the calorimeter (fake 
Kl -crash), while the Kl meson decays close enough to the interaction point to produce six 
prompt photons [5j. To suppress this kind of background we first reject events with charged 
particles coming from the vicinity of the interaction region. Next, taking advantage of the 
differences in the /?* and E cr distributions between the tagged Kg events and the fake Kl~ 
crash, we have tightened cuts on these variables: E cr > 150 MeV and 0.20 < /?* < 0.225. 
This improves by a factor 12 the rejection of this background with respect to the previous 
analysis [3] . The second source of background originates from wrongly reconstructed K$ — > 2-7T 
decays. The four photons from this decay can be reconstructed as six due to fragmentation 
of the electromagnetic showers (splitting). These events are characterized by one or two low- 
energy clusters reconstructed very close to the position of the genuine photon interaction in the 
calorimeter. Additional clusters can come from the accidental time coincidence between the (ft 
decay and machine background photons from DA<I>NE [S]. In the next stage of the analysis 
we select only kinematically well defined events. To this end we perform the kinematical fit 
procedure based on the least squares method imposing energy and momentum conservation, the 
kaon mass and the velocity of the six photons in the final state. The x 2 distribution of the fit for 
data and background simulation, Xfui ls shown in Fig. together with the expected distribution 
for signal events. Cutting on Xju considerably reduces the background from bad quality 
reconstructed events while keeping large signal efficiency [5]. In order to reject events with split 
and accidental clusters we look at the correlation between two x 2 ~hke discriminating variables 
C2tt and C3?r- (2n is calculated by an algorithm selecting four out of six clusters best satisfying 
the kinematic constraints of the two-body decay, therefore it verifies the K$ — > 27r° — > 47 
hypothesis. The pairing of clusters is based on the requirement m 77 = m^o , and on the opening 
angle of the reconstructed pions trajectories in the K$ center of mass frame. Moreover, we check 
the consistency of the energy and momentum conservation in the (ft ~^ KsKl, K$ —> 27r° decay 
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Figure 3. Distribution of Xfu f° r 
the tagged six-photon sample for data 
( black points), background simulation 
(solid histogram), and simulated K$ — > 
3tt° signal (dashed histogram). 
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Figure 4. Distributions of events in the C37r-C27r plane, for six-photon sample tagged by K^— 
crash for data (left), and for the simulated Kg — > 3ir° decays (right). The solid lines show the 
signal region definition chosen using the optimization procedure described in the text. 



hypothesis. The X3V instead verifies the signal hypothesis by looking at the reconstructed masses 
of three pions. For every choice of cluster pairs we calculate the quadratic sum of residuals 
between the nominal it mass and the invariant masses of three photon pairs. As the best 
combination of cluster pairs, we take the configuration minimizing £3^. The distributions in the 
C37T— plane for the data and K$ — > 37r° simulated signal are shown in Fig. 01 Signal events are 
characterized by small values of £3^ and relatively high fan [H]> which indicates the choice of the 
signal region definition (signal box). In order to improve the quality of the photon selection using 
we cu t on the variable AE = (m$/2 — J2 Ey t )/aE where 7« stands for the i-th photon from 
four chosen in the xItt estimator and oe is the appropriate resolution. For K$ — > 2tt° decays plus 
two background clusters, we expect AE ~ 0, while for K$ — > 3tt° AE ~ tti^o/ije- To further 
reject surviving K$ — > 2ir° events with split clusters, we cut on the minimal distance between 
centroids of reconstructed clusters, R m in, considering that the distance between split clusters is 
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Figure 5. Distributions of AE and R m i n discriminating variables for six-photon events for data 
(black points) and background simulations (red histograms). The dashed histograms represents 
simulated K$ — > 37r° events. 
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Figure 6. Distributions of Xfu f° r six-photon events in the signal box (left) and AE for 
six-photon events in the signal box applying the x}u < ^7.2 cut (right). Black points are 
data, background simulation is the red histogram. The dashed histogram represents simulated 
Ks — > 37T° events. 



on average smaller than the distance between clusters originating from 7's of Ks —> 3tt decay. 
Distributions of these two discriminant variables are presented in Fig. [5j Before counting of 
the Ks —> 37T° candidate events, we have optimized the event selection in order to reduce the 
background as strongly as possible while keeping high signal efficiency. Cuts on the discriminant 
variables have been refined minimizing f C ut{Xfu> C2w, (3^, A.E, Rmi n ) = N up /e^, where £371- stands 
for the signal selection efficiency and N up is the mean upper limit (at 90% C.L.) on the expected 
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Figure 7. Distributions of Rmin for six-photon events in the signal box applying the Xju < 57.2 
cut (left), and applying xju < 57.2 and AE > 1.88 cuts (right). Black points are data, 



background simulation is the red histogram. 
Kg — >• 37T° events. 



The dashed histogram represents simulated 



number of signal events. It is calculated on the basis of the corresponding expected number 
of background events B exp = B exp (xf it ,C2n,(3n,AE,R m i n ) from simulation pp. As the result 
of the optimization we have obtained the following values for cuts on discriminant variables: 
Xfu < 57.2, AE > 1.88 and R m i n > 65 cm. The signal box is defined as: 4 < Qi-k < 84.9 
and C3tt < 5.2 |5j. Since the expected number of background events was estimated using the 
Monte Carlo simulations, we have checked at each stage of the analysis how well the simulation 
describes the experimental data. Distributions of AE and Rmin variables are presented in 
Fig. [6] and Fig. [7] for events in the signal box. In the right panel of Fig. [7] we present also the 
Rmin distribution just before the last cut R m in > 65 cm marked with a black arrow. One can 
see a good agreement between simulations and measured data at every stage of the analysis. At 
the end of the analysis we have found zero candidates in data and in the simulated background 
sample. 

As it was mentioned, the Ks — > 27r° normalization sample is selected requiring four prompt 
photons. The Monte Carlo simulation shows a negligible (0.1%) amount of (j) — > K + K~ 
background events, thus no further discriminant analysis was required. In the same tagged 
sample of Ks mesons we have found iV^ = (7.533 ± 0.018) • 10 7 events. Using simulations we 
have also determined the Ks — > 2tt° — > 47 selection efficiency: 62^ = 0.660 ±0.002 S £ a £. The final 
number of produced Ks — > 27r° events amounts to: N norra = N2 7T /e2-K = (1.142 ± 0.005) • 10 8 . 

4. Results 

With cuts defined in the previous section at the end of the analysis we count candidates with 
background events expected from Monte Carlo simulated with an effective statistics of two times 
that of the data. Hence, we have estimated an upper limit on the Ks — > 3-7T branching ratio 
taking into account systematic uncertainties related to the number of background events and to 
the determination of the acceptance and selection efficiencies for the signal and normalization 
samples. The detailed description of the evaluation of systematic uncertainties can be found in 
Ref . [U [5] , here we would like only to stress that they are less than 5% and are negligible in the 



calculation of the limit. 

In the conservative assumption of no background, we estimate an upper limit on the K$ — >• 37T° 
branching ratio at 90% confidence level: 

BR{K S -»• 3vr°) < 2.6 • 10~ 8 , (1) 

which corresponds to an improvement by a factor of about 5 compared to our previous search [S]. 
This result can be translated into a limit on 1 77000 1 : 

1 77000 1 = \A(Ks -» 3tt°)/A(K l ->■ 37r°)| < 0.0088 at 90% C.L. (2) 
5. Summary and outlook 

As a result of the KLOE data set analysis, gathered in the 2004-2005 data taking period, we 
have set the new upper limit for the K$ — > 37r° branching ratio at the 90% confidence level, 
which is almost five times lower than the latest published result [3j. However, the search for 
the K$ — > 37T° decay will be continued by the KLOE-2 collaboration, which is continuing and 
extending the physics program of its predecessor [11]. For the forthcoming run the KLOE 
performance have been improved by adding new subdetector systems: the tagger system for the 
77 physics, the Inner Tracker based on the Cylindrical GEM technology and two calorimeters 
in the final focusing region [12]. These new calorimeters will increase the acceptance of the 
detector, while the new inner detector for the determination of the K$ vertex will significantly 
reduce the contribution of the background processes involving charged particles. Increasing the 
statistics and acceptance of the detector while significantly reducing the background gives a 
realistic chance to observe the K$ — > 37r° decay for the first time in the near future. 
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